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Abstract. We consider manifolds with conic singularites that are iso- 
metric to 1" outside a compact set. Under natural geometric assump- 
tions on the cone points, we prove the existence of a logarithmic resonance- 
free region for the cut-off resolvent. The estimate also applies to the 
exterior domains of non-trapping polygons via a doubling process. 

The proof of the resolvent estimate relies on the propagation of sin- 
gularities theorems of Melrose and the second author [23] to establish a 
"very weak" Huygens' principle, which may be of independent interest. 

As applications of the estimate, we obtain a exponential local energy 
decay and a resonance wave expansion in odd dimensions, as well as a 
lossless local smoothing estimate for the Schrodinger equation. 



1. Introduction 

In this paper we consider a manifold X of dimension n with conic singu- 
larities that is isometric to M. n outside a compact set. We impose geometric 



hypotheses (elucidated in Section 2.4 as Assumptions 1-3) that 

(1) The flow along "geometric" geodesies is non-trapping. (Geomet- 
ric geodesies are those that miss the cone points or that are every- 
where locally given by limits of families of geodesies missing the cone 
points.) 

(2) No three cone points are collinear. 

(3) No two cone points are conjugate to each other. 

Our main result is as follows (throughout the paper, A denotes the Laplacian 
with positive spectrum): 

Theorem 1. For \ e C^f(X), there exists 5 > such that the cut-off resol- 
vent 

X(A-A 2 )- 1 X 

can be analytically continued from Im A > to the region 
ImA > -5 log Re A, Re A > <5 _1 
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and for some C, T > enjoys the estimate 
in this region. 



X (A-A 2 )- 1 xL,w,^G|Ar 1 e r|ImA| 



As shown by Lax-Phillips [20] and Vainberg [29] in certain geometric set- 
tings and later generalized by Tang-Zworski ]28] to "black-box" perturba- 
tions, if the dimension n is odd, then Theorem [T] results in a decay estimate 
for solutions to the wave equation in such a geometry, and indeed in a full 
resonance-wave expansion for solutions to the wave equation. Let T> s denote 
the domain of A s l 2 (see Section [2] below) and let sin ty/~A/y/A~ be the wave 
propagator. Let x equal 1 on the set where X is not isometric to M. n . 

Corollary 1. Let n be odd. For all A > 0, small e > 0, t > sufficiently 
large, and f e T>\, 



* A f eRes(A) m=0 

Im A>-A 

where the sum is of resonances of A, i.e. over the poles of the meromorphic 
continuation of the resolvent, and the Wj tTn are the associated resonant states 
corresponding to Xj. The error satisfies 

\\E A (t)\\ v ^ L2 ^C e e-^. 

In particular, since the resonances have imaginary part bounded above by 
a negative constant, x sin ^^ xf ^ s exponentially decaying. 

(We refer the reader to Theorem 1 of [28j for details of the resonance wave 
expansion.) 

Another consequence of our resolvent estimate is a local smoothing es- 
timate without loss for the Schrodinger equation. Local smoothing esti- 
mates were originally established for the Schrodinger equation on K n by 
Sjolin [23], Vega [HU], Constantin-Saut [TT], Kato-Yajima [18] . and Ya- 
jima [32] . Doi [15] showed that on smooth manifolds the absence of trapped 
geodesies is necessary for the local smoothing estimate to hold without loss. 
We show that even in the presence of very weak trapping due to the diffrac- 
tive geodesies, the local smoothing estimate holds without loss. 

Corollary 2. Suppose u satisfies the Schrodinger equation on X: 

r l d t u(t, z) + Au(t, z) = 

u(0,z) = u {z) e L 2 {X) 

Then for all x e C^P(X), u satisfies the local smoothing estimate without 
loss: 



T 

X«(*)llx>i/2 dt sS C T \\uo\\ L 2( X ) 
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This result follows directly from our Theorem[T]by an argument of Burq [3]. 

Another application of the resolvent estimate of Theorem [T] is to the 
damped wave equation. Although we do not pursue it here, under suitable 
convexity assumptions (e.g., if no geodesic passing through the perturbed 
region re-enters it; see Datchev-Vasy \12\ 113] for more general conditions), 
it is possible to obtain decay estimates for the damped wave equation on 
conic manifolds when the only undamped geodesies are diffractive ones. 
This relies on a gluing construction of Datchev-Vasy to obtain a suitable 
resolvent estimate and on the recent work of Christianson, Schenck, Vasy 
and the second author [9] to yield the estimate. 

In addition to applying to manifolds with cone points, our results also ap- 
ply to the more elementary setting of certain exterior domains to polygons 
in the plane. Let c K 2 be a compact region with piecewise linear bound- 
ary. We further suppose that the complement M 2 \f2 is connected, that no 
three vertices of Q are collinear, and that IR 2 \Q is non-trapping, in the sense 
that all billiard trajectories not passing through the vertices of $7 escape to 
infinity^] Figure [l] illustrates an example of such an exterior domain. For 
this class of domains, the analogue of Theorem [T] holds. 




Figure 1 . An example of a domain to which Corollary [3] 
applies. The dashed line represents a trapped diffractive or- 
bit. 



Corollary 3. If X = M 2 \f2 is the exterior of a non-trapping polygon with 
no three vertices collinear and A is the Dirichlet or Neumann extension of 
the Laplacian on X, then the result of Theorem^ holds for the resolvent on 
X. 

The proof of Corollary [3] relies on reducing the problem to one on a surface 
with conic singularites. Indeed, such an exterior domain can be doubled by 
gluing together two copies of it across the common boundary; this results 
in a manifold with cone points, corresponding to the vertices of the initial 
polygonal domain, and with two ends each isometric to M 2 . Solutions to the 
wave equation in this "doubled" manifold are closely related to solutions to 
the Dirichlet or Neumann problem on the original exterior domain via the 

^In fact we require a slightly different condition. We ask that all billiard trajectories 
that are locally approximable by trajectories missing the vertices escape to infinity. This 
is not quite the same condition but is generically equivalent. 
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method of images. Our results hold for the exterior problem to such non- 
trapping polygons as well, although this entails some mild complication in 
the proof (the introduction of "black-box" methods) — see Section [5] below. 
In particular, our result affirmatively answers a conjecture of Chandler- 
Wilde, Graham, Langdon, and Spence [3j. In the case when the obstacle is 
star-shaped, we remark that the exponential energy decay is a consequence 
of the classical technique of Morawetz estimates (see e.g. Lemma 3.5 of 
[5])); we believe that the estimate for general non-trapping polygons is new, 
however. 

We note that stronger estimates than those of Theorem [lj are known to 
hold in the case of a non-trapping metric or even an appropriately non- 
trapping "black box" perturbation such as a smooth non-trapping obstacle 
(see |20j . [29], |28j): in these cases there are finitely many resonances above 
any logarithmic curve ImA > —N log Re A. That the result here is likely 
to be sharp can be seen from the explicit computation of Burq [2], who 
shows that in the case of obstacle scattering by two strictly convex analytic 
obstacles in M 2 , one of which has a corner, the resonances (poles of the 
analytic continuation of the resolvent (A — A 2 ) -1 ) are located along curves 
of the form ImA = —Clog Re A. Burq's setting is not of course exactly that 
of manifolds with cone points, but is suggestively close to that of polygonal 
domains discussed above. (Similar logarithmic strings of resonance poles 
also appear in Zworski |33j where they are shown to arise from finite order 
singularities of a one-dimensional potential, substantiating heuristics from 



By contrast, what seems the weakest trapping possible in the setting of 
smooth manifolds, a single closed hyperbolic geodesic, is known in certain 
settings to produce strings of resonances along lines of constant imaginary 
part [TO] , and hence yields an analytic continuation to a smaller region than 
that shown here, which does not permit a resonance wave expansion in the 
strong sense of Corollary [I] except in very special cases [8] . 



The fact that the estimates demonstrated here are weaker, by only a 
very small margin, than those for non-trapping situations, reflects that fact 
that cone points induce a kind of "weak trapping:" there exist geodesies 
connecting every pair of cone points, and concatenation of such geodesies 
starting and ending at the same cone point should be considered a legiti- 
mate geodesic curve in a conic geometry. In particular, such concatenations 
of geodesies are known (generically) to propagate singularities of the wave 
equation on exact cones by results of Cheeger- Taylor [61 [7J; this reflects 
the diffraction of singularities by the cone point. Melrose and the second 
author |23| subsequently showed that on any manifold with conic singular- 
ities, the propagation of singularities is limited to geodesies entering and 
leaving a given cone point at the same time ( "diffractive propagation"). It 
was further shown in [23j that the fundamental solution of the wave group 
with initial pole near a cone point was smoother along generic geodesies 
emerging from the cone point than along those that are approximable by 
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geodesies emanating from the initial pole and missing the cone point; this 
"smoothing effect" in fact holds for any solution that satisfies an appropriate 
nonfocusing condition with respect to the cone point in question (see Sec- 



tion 2.2 below). Thus, colloquially, [23j showed that "diffracted singularities 
are smoother than geometrically propagated singularities." It also showed 
that the spherical wavefront of diffracted singularities is a conormal wave. It 
is the smoothing property and the conormality that play an essential role in 
the proof of Theorem [TJ The proof proceeds via another result which may be 
of independent interest, a theorem on the weak non-trapping of singularities 
for manifolds with cone points. In the following theorem, U(t) denotes the 
wave group, and £ r denotes the Sobolev space of energy data V r ©P r _i. 

Theorem 2. Let % £ For any set, there exists T s » such that 

whenever t > T s , 

xU(t)x : £ r -> £r+s 

for all r. 

We recall that Huygens' Principle, valid in odd dimensional Euclidean 
space, says that xU(t)x is eventually identically zero. More generally, in 
even dimensional Euclidean space or indeed in any "non-trapping" metric 
in which all geodesies escape to infinity, xU{t)x eventually has a Schwartz 
kernel in C 00 . Our Theorem [2] is weaker yet: here the cut-off wave kernel is 
as smooth as one likes, after a sufficiently long time. 

Resolvent estimates similar to ours have been previously demonstrated by 
Duyckaerts [16] for operators of the form A + V where A is the Euclidean 
Laplacian and V has multiple inverse-square singularities: these singularities 
are analytically similar to (albeit geometrically simpler than) cone points. 

2. Geometric set-up 

The basic material in this section on conic geometry comes from [23J 
while the more detailed discussion of the global geometry of geodesies is 
taken from |31j . 

Let X be a noncompact manifold with boundary, K a compact subset of 
X, and let g be a Riemannian metric on X° such that X\K is isometric 
to the exterior of a Euclidean ball M. n \B n (0, Rq) and such that g has conic 
singularities at the boundary of X: 

g = dx 2 + x 2 h(x, dx, y, dy); 

here g is assumed to be nondegenerate over X° and h\gx induces a metric 
on dX. We let Y a , a = 1,...N denote the components of dX; we will 
refer to these components in what follows as cone points, as each boundary 
component is a single point when viewed in terms of metric geometry. 
We further let 

M = Rx X 
denote our space-time manifold. 
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We recall from Theorem 1.2 of [23] that by judicious choice of coordinates 
x,y on a collar neighborhood of dX, we may reduce g to the normal form 

(1) g = dx 2 + x 2 h(x,y,dy), 

where h is now a family (in x) of metrics on Y. Then the curves y = const, are 
geodesies, with x the length parameter. Indeed, the curves of this form are 
the only geodesies reaching dX and they foliate a neighborhood of dX. We 
let x J-& denote the collection of the continuations of forward and backward 
bicharacteristics in T*X° which reach the boundary component Y a in time 
\t\ ^ s (with T denoting "flowout" of the cone point Y a ). Thus for small s, 
in canonical coordinates £, rj dual to x, y, 

x F s a = U n {x e {0,s),y e Y a ,£ e R, v = 0} 

where U is a neighborhood of the single boundary component Y a containing 
a component of x < s. We further refer to points in X T^ as incoming or 
outgoing with respect to the cone point according to whether they reach 
the boundary at positive or negative time respectively under the flow (this 
separates the manifold into components). We will also be concerned with the 
corresponding flowout sets in space-time. Letting £ denote the characteristic 
set of □ = D\ - A on T*M° we define 

K = {(t, r, z, e £ : (z, e X K) <= T* M° . 

As discussed in [22] (where the notation T was first used) the manifolds 
X F&, J~a are coisotropic conic submanifolds of T*X°, resp. T*M°. 

We let T> s denote the domain of the s/2 power of the Friedrichs extension 
of the Laplacian on Cf{X°). Note that this agrees with the ordinary Sobolev 
space H s away from the cone points (and was characterized in [231 Section 
3] in terms of the scale of weighted b-Sobolev spaces). Let 

S s = V S ®V S - X 

denote the corresponding space of Cauchy data for the wave equation, and 
let 

U(t) = expit 
denote the wave propagator, hence 

U(t) :£ s ^£s 

for each s e K. We will frequently need to deal with error terms that are 
residual in the scale of space £ s , so we define 

with the additional c subscript denoting compact support in X. In dealing 
with wave equation solutions as functions in spacetime, it is convenient to 
think of them lying in the Hilbert space 

L 2 ([0,f];£ s ) 



(I o 
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Figure 2. A family of geodesies missing the cone point (i.e., 
boundary component), limiting to a pair of geodesies entering 
and leaving it normally together with a geodesic of length it 
in the boundary connecting the two. 

with T » T s taken large enough to encompass all time intervals under con- 
sideration. We thus denote this space 

L 2 £ s 

for brevity. We recall that solutions to the wave equation in L 2 £ s have 
unique restrictions to fixed-time data lying in £ s , and will use this fact 
freely in what follows. 

For convenience, we will equip the cosphere bundle S*X with a Riemann- 
ian metric inducing a distance function, denoted 

2.1. Geometric and diffractive geodesies. We now recall the different 
notions of "geometric" and "diffractive" bicharacteristic which enter into 
the propagation of singularities on manifolds with cone points. 

Definition 1. A diffractive geodesic on X is a union of a finite number of 
closed, oriented geodesic segments 71, ... , in X such that all end points 
except possibly the initial point in 71 and the final point of 7 at lie in Y = dX , 
and 7i ends at the same boundary component at which begins, for 
i = 1,...,JV-1. 

A geometric geodesic is a diffractive geodesic such that in addition, the 
final point of 7^ and the initial point of 7^+1 are connected by a geodesic 
of length 7r in a boundary component Y a (w.r.t. the metric ho = h\y a ) for 
i = 1,...,JV-1. 

The proof the following proposition was sketched in |23j, and yields the 
equivalence of the above definition of "geometric geodesic" with the more 
casual one used in the introduction above: 

Proposition 1. The geometric geodesies are those that are locally realizable 
as limits of families of geodesic in X° as they approach a given boundary 
component. 

Remark 1. We note that while every geometric geodesic is locally approx- 
imable by smooth geodesies in X°, a geodesic undergoing multiple inter- 
actions with cone points may not be globally approximable in this sense. 
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Figure 3. A non-approximable geometric geodesic in the 
triply-slit-plane (or its branched cover), with an approximat- 
ing broken geodesic in X°. 

Figure [3] shows such a situation. It is most easily interpreted as showing 
a domain with boundary given by the three slits; geodesies then reflect off 
smooth parts of the boundary, and geometric geodesies either pass straight 
through the end of the slits or reflect specularly as if the slits continued. 
Then the vertical line, while it can be uniformly approximated by broken 
geodesies as shown, cannot be approximated globally: any approximating 
geodesic would have to reflect off one or another of the slits. To place this 
example in the context of the current paper rather than that of domains 
with boundary, we should instead interpret the picture as showing one sheet 
of a two-sheeted ramified cover of M 2 , with the slits representing branch 
cuts. This makes the ends of the slits into cone points, with the link of each 
cone point a circle of circumference 47T. In this situation, any unbroken ap- 
proximating geodesic would have to move onto the other sheet of the cover 
by passing through one of the slits, hence could not globally approximate 
the line shown, which remains on a single sheet of the cover. 

This very simple example has three collinear cone points, which we are 
ruling out by hypothesis; however one can make other examples involving 
only two interactions with cone points, and this is permitted by our geomet- 
ric hypotheses. For this reason Assumption [T] is formulated so as to cover 
geometric geodesies explicitly, rather than just as a uniform statement on 
geodesies in X° (which might be preferable). Ruling out propagation along 
non-approximable geometric geodesies, if indeed true, would likely require 
somewhat delicate second-microlocal arguments. 

As we will need to understand the lift of the geodesic flow to the cotangent 
bundle, it is helpful to see how this can be accomplished uniformly up to d X 
(although in this paper, microlocal considerations will only arise over X°, 
which is a considerable simplification). We let b T*X denote the b-cotangent 
bundle of X, i.e. the dual of the bundle whose sections are smooth vector 
fields tangent to dX. Let b S*X denote the corresponding sphere bundle. Let 
£ dx/x + 7]-dy denote the canonical one-form on h T* X. (We refer the reader 
to Chapter 2 of |21j for a further explanation of "b-geometry," of which we 
only use the rudiments here.) 
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Let Kg be the Hamilton vector field (with respect to the symplectic form 
d(£dx/x + rj ■ dy)) for g/2 = (£ 2 + h(x, y, r]))/(2x 2 ), the symbol of A/2 on 
b T*X; note that K g is merely the geodesic spray in b T*X with velocity 
*Jg. It is convenient to rescale this vector field so that it is both tangent to 
the boundary of X and homogeneous of degree zero in the fibers. Near a 
boundary component Y a , for a metric in the reduced form Q, we have (see 



H Ya {x) + (e 2 + h(x,y, V ) + -—)d ii + ixd x V 

where Hy a {x) is the geodesic spray in Y a with respect to the family of 
metrics h(x,-). Hence the desired rescaling is 

Z = ^-Kg. 

(Note here that g refers to the metric function on the cotangent bundle 
and not the determinant of the metric tensor.) By the homogeneity of 
if we radially compactify the fibers of the cotangent bundle and identify 
b S*X with the "sphere at infinity" then Z is tangent to h S*X, and may 
be restricted to it. Henceforth, then, we let Z denote the restriction of 
(x/*Jg)K g to the compact manifold b S*X on which the coordinates £, 77 
have been replaced by the (redundant) coordinates 



Z vanishes only at certain points x = fj = over dX, hence the closures of 
maximally extended integral curves of this vector field can only begin and 
end over dX. Since Z is tangent to the boundary, such integral curves either 
lie entirely over dX or lie over dX only at their limit points. Interior and 
boundary integral curves can meet only at limit points in {x = fj = 0} a 
b S*X. 

It is helpful in studying the integral curves of Z to introduce the following 
way of measuring their lengths: Let 7 be an integral curve of Z over X° . 
Let k denote a Riemannian metric on b S*X° such that k(Z, Z) = 1. Let 

(3) u = xk{-,Z) en l { b S*X). 

Then 

a; = xk(dj/ds, Z) ds = ——k(K g ,Z)ds= ds = length (7) 

where s parametrizes 7 as an integral curve of Kg/^/g, the unit speed geo- 
desic flow. With this motivation in mind, we now define, for each t e M.+ , 
two relations in b S*X, a "geometric" and a "diffractive" relation. These 
correspond to the two different possibilities for geodesic flow through the 
boundary. 
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Definition 2. Let p,ge b S*X. We write 

G,t 
p~q 

if there exists a continuous, piecewise smooth curve 7 : [0, 1] — b S*X with 
7(0) = p, 7(1) = q, such that [0, 1] can be decomposed into a finite union of 
closed subintervals Ij, intersecting at their endpoints, where 

(1) on each I?, 7 is a (reparametrized) positively oriented integral curve 
of Z in b S*X, 

(2) On successive intervals Ij and Ij+i, interior and boundary curves 
alternate, 

(3) ^oj = t, with ui as defined in ([3]). 

We write 

D,t 
p~q 

if there exists a piecewise smooth (not necessarily continuous) curve 7 : 
[0, 1] — b S*X with 7(0) = p, 7(1) = q, such that [0, 1] can be decomposed 
into a finite union of closed subintervals Ij, intersecting at their endpoints, 
where 

(1) on each J?, 7 is a (reparametrized) positively oriented integral curve 
of Z in b S*X°, 

(2) the final point of 7 on Ij and the initial point of 7 on Ij+i lie over 
the same component Y a of dX, 

(3) S 7 ^ = t. 

Integral curves of Z over X° are lifts of geodesies in X°, and it follows 
from Q that the maximally extended integral curves of Z in b Sg X X are 
lifts of geodesies of length ir in dX (see [23] for details), hence: 

G t 

Proposition 2. p~q iff p and q are connected by a (lifted) geometric geo- 
desic of length t. 

p^~q iff p and q are connected by a (lifted) diffractive geodesic of length 

t. 

A very important feature of these equivalence relations, proved in [31], is 
the following: 

Q £ 

Proposition 3. ([31]. Prop. 4) The sets {(p,q,t) : p~<?} and {(p,q,i) : 
p D ~q) are closed subsets of b S*X x b S*X x R+. 

2.2. Propagation and diffraction of singularities. We now recall the 
key propagation results from [23 . In order to do this, we must introduce 
the notions of coisotropic regularity and nonfocusing with respect to a conic 
coisotropic submanifold X of T*M°. We will usually take X = T a ~ J-^ with 
S taken larger than all times on which we will be considering propagation 
of singularities (and hence ignored in our notation). 
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Let A4 denote the space of compactly supported pseudodifferential oper- 
ators Ae ^ ,1 (M°) that are characteristic on X. Let A be the filtered algebra 
generated by M. over *°(M), with A k = A n ^ k {M). Fix a Sobolev space 
U. Fix a compact set K cz S*M°. 

We now recall the following definitions from |22j. 

• We say that u has coisotropic regularity with respect to X of order k 
relative to ~H in K if there exists ^4 e fy°(M), elliptic on K and sup- 
ported over M° such that A k Au c We say that u has coisotropic 
regularity relative to % in K if it has coisotropic regularity of order 
k for all kl 

• We say that u satisfies the nonfocusing condition of degree k with 
respect to X relative to H on K if there exists A e $°(M), elliptic 
on if and supported over M° such that cz „4 fe %. We say that u 
satisfies the nonfocusing condition relative to T~L on K if it satisfies 
the condition to some degree. 

If K is omitted (which will only be the case when the distribution is micro- 
supported away from the boundary), the relevant condition is assumed to 
hold on all of S*M° . 

In the special case when X = T a and we work in a collar neighborhood 
of the boundary component Y a in which the metric has the form ([!]), these 
definitions simplify considerably, as the module Ai is then generated by the 
operators d Vj . A distribution u microsupported over such a neighborhood is 

coisotropic of order k relative to H iff (Ay) fe//2 u e H, while it satisfies the 
nonfocusing condition relative to T~L iff there exists N with (Ay) - u e %. 
(Here Ay denotes the Laplacian in the y- variables with respect to the family 
of metrics h(x) on Y a .) 

We are now in a position to recall the main results of |23j . (Our notation 
however sticks more closely to that of [22], which treats the more general 
case of edge manifolds.) 

Note that here and henceforth we employ the notation s — to mean s — e 
for every e > 0. 

Proposition 4. Let u e C(R;D r ) be a solution to the wave equation on M. 
Fix a point p e T a outgoing with respect to the cone point Y a . 

(1) If u is microlocally in H s on all incoming bicharacteristics in J- a 
that are diffractively related to p then u e H s microlocally at p. 

(2) Assume that u is nonfocusing with respect to H s on J- a and has no 
wavefront set along incoming bicharacteristics in X a that are geo- 
metrically related to p. Then u e H s ~° microlocally at p, and enjoys 
coisotropic regularity relative to H s ~° in a neighborhood of p. 

Remark 2. 

• The coisotropic regularity part of this result is left slightly implicit 
in the results of [23] but follows easily by interpolating part (iii) of 
Theorem 8.1 of that paper, which gives coisotropic regularity relative 
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to some Sobolev space, with the overall regularity of the solution 
microlocally near p. 
• The following consequence is more germane to what follows (and 
perhaps easier to digest): Fix s, k with r < k < s. Say our solution 
is nonfocusing with respect to H s . If a singularity in WF fc (the set 
where u $ H k microlocally) arrives at Y a along just one ray in J- a , 
with u microlocally smooth along all other arriving rays, then it 
may produce milder singularities, at worst in H s ~° and coisotropic 
relative to this space, along rays diffractively related to the incoming 
singularity, but has strong singularities in WF fc along (at least some 
of) those geometrically related to it. 

The crucial example is of course the fundamental solution 

sintVA 



with p close to a boundary component Y a . This solution is overall in jj- n / 2 + 1_ 
but it satisfies the nonfocusing condition relative to H 1 ^ 2 ^ . When the singu- 
larity strikes Y a is produces strong singularities (no better than }j^ n l 2 + l - Q ^ 
along the geometric continuations of the geodesic from p to Y a but only 
weaker singularities in ff 1 / 2-0 along the rest of the points in a spherical 
wavefront emanating from Y a , along which the solution in fact enjoys La- 
grangian regularity. 

The same reasoning that applies to the fundamental solution in fact ap- 
plies to a solution given by a Lagrangian distribution with respect to a 
Lagrangian manifold intersecting J- a transversely. In this paper, however, 
we will be concerned with a slightly different setting, in which our hypothe- 
ses on the solution are themselves that of coisotropic regularity, but with 
respect to a different coisotropic manifold, intersecting T a transversely (in 
particular, with respect to for some (3). In applying Proposition [4] above, 
we will need a result in order to verify the nonfocusing hypotheses. This 
result is discussed in the next section. 



2.3. Intersecting coisotropics. In this section we prove the following re- 
sult allowing us to verify the nonfocusing hypotheses of Proposition [4} 

Proposition 5. Let X denote a conic coisotropic manifold of codimension 

n — 1 intersecting J- a in such a way that the isotropic foliations of the two 

coisotropic manifolds are transverse. Let u enjoy coisotropic regularity with 

respect to I, relative to H s . Then u is nonfocusing with respect to J- a relative 
to H s+(n-l)/2-a_ 

The proof of the proposition relies on two lemmas. We first show it in a 
model setting: 
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Lemma 1. In R n i+ n 2 with coordinates z e R ni , z' e W 12 , let X, X denote 
the two model coisotropic manifolds 

1 = {C = 0}, X = {z = 0}. 

If a compactly supported distribution enjoys coisotropic regularity with re- 
spect to X relative to H s then it is nonfocusing with respect to X relative to 
£P+ni/2-o_ That is to say, there exists N » such that 

(A z y N u e Hl + c ni/2 -° 

Note that n\ is the dimension of the leaves of the isotropic foliation. 

Proof. By applying powers of A we may reduce to the case s = 0. Coisotropic 
regularity with respect to X relative to L? means that 

z a D^,Dju e L 2 , if |a| = \p\ + | 7 |, 

as the vector fields ZiD Zj and ZiD z ^ have symbols cutting out X. 

In particular, the iterated regularity under ZiD Zk means that we have 

z a Z)>eL 2 (M™ 2 ;/ (0) (M ni ;0)), |a| = 

where I'M (K m ; 0) denotes the space of functions v e H s (R ni ) enjoying it- 
erated regularity under vector fields ZiD Zk , i.e., conormal regularity with 
respect to the origin. 

Now we claim that if p < n\/2 is positive, 

\ Z \-PI^ {R ni ;0) c H- ni / 2 -°{R ni ). 

Indeed, interpolation shows that for a compactly supported u e /( ) (R ni ; 0) , 

Af \z\ s ueL 2 

for all s > 0. Hence Sobolev embedding shows that 1^ is contained in 
|^|-m/2-0 L co_ ThuS; for p < ni y 2) |^|-Pj(0)( R n 1;0 ) ig contained in \ z \- n ^L x 

for some e > 0. This in turn implies that for compactly supported u e 
|zp p /(°)(M ni ;0), u is bounded and so <0 _ni/2 ~ e n e £ 2 for all e. Taking the 
inverse Fourier transform yields u e _£f~ ni//2 ~ <E (IR ri1 ) for all e > 0. 

Thus, applying powers of z a Dj,u (and again interpolating to deal with 
fractional powers) we find that 

ue H s (R n2 ;H- ni / 2 -°(R ni )), s < m/2. 

Finally, this implies that for N > ni/2, 

(A z y N ue H s (R n2 ;H n ^ 2 {R ni )), s < m/2. 

Since H ni /2-0( R n 2 . iW 2 (IR ni )) c F ni / 2 -°(K ni+n2 ) this implies nonfocus- 
ing with respect to n\/2 — as desired. □ 

Lemma [T] is in fact quite general, as the form of the coisotropic distri- 
butions employed there is a normal form for intersecting conic coisotropic 
manifolds with transverse foliations: 
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Lemma 2. Let X,X be conic coisotropic submanifolds o/T*(R n ), each of 
codimension k < n, and intersecting at p e T*(R n ) in such a way that the 
isotropic foliation of each is transverse to the other. Then there exist local 
symplectic coordinates (z, z' , C> C) * n which p lies at the origin, and 

1= {C = 0}, X = {z = 0}. 

Proof. We choose a degree-one homogeneous function / vanishing simply 
along X. Thus its Hamilton vector field H f has degree zero and is by definition 
tangent to X and transverse to X. On the flowout of X by H f we let g denote 
the homogeneous, degree-zero function measuring the time of flowout from 
X; we may further extend g to a homogeneous degree-zero function (also 
denoted g) on T*M n such that 

{/,$} = ! 

and with 

g = 0, dg =£ on X 

(To make this extension, take coordinates, not necessarily symplectic, in 
which X lies along coordinate axes, and define g locally to be the time to 
flow out to a nearby point from a hyperplane containing X and transverse 
to H/.) 

Now by Theorem 21.1.9 of |19j . we may extend the coordinates 

{zkXk) = (g,f) 

to a full system of homogeneous symplectic coordinates. Thus, we have 
locally achieved 

Xc{C fe = 0}, l^{z k = 0}. 

Since H^ fc ,H 2fc are respectively tangent to X,X, these manifolds are locally 
products of the form 

1' x {{z k e R,( k = 0)}, X x {{z k = 0,C fe 6 M)} 

with X',X' coisotropic submanifolds of T*IR n_1 satisfying the hypotheses of 
the lemma with n replaced by n — 1 and k replaced by k — 1. 

The result then follows by induction. □ 

Proof of Proposition [5]. Proposition[5]now follows from the two lemmas above: 
we can bring the two coisotropic manifolds J- a and X into the normal form 
given by Lemma [2j Lemma [T] then shows that nonfocusing holds, with 
n\ = n — 1, the dimension of the isotropic foliation (which in the case of J- a 
has leaves obtained by varying y e Y a ). □ 
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2.4. Geometric Assumptions. We are now able to state our geometric 
assumptions: 

Assumption 1. Let VL 3 K be an open set with X\VL ~ M. n \B n (0, R{) 
for some R± » 0. We assume that there exists a time Tq > such that any 
geometric geodesic starting in S^X° leaves f2 in time less than Tq ( and does 
not come back, by convexity of the ball). 

Assumption 2. No geometric geodesic passes through three cone points. 

Assumption 3. No two cone points Y ai Yp are conjugate to one another 
along geodesies in S*X° of lengths less than Tq in the sense that whenever 
s + t < To, J~a and J-p intersect transversely for each a, f3. 

Remark 3. Assumption [I] is a quantitative statement of non-trapping of 
geodesies that do not hit the cone points. 

Assumption [2] is generic at the formal level of dimension counting, as a 
geodesic arriving at Y~ from Y a can be geometrically continued from a family 
of points in Y^ of dimension n — 2 (those points at distance ir in Y^ from 
the arrival point); on the other hand there is a (generically) discrete set of 
departure points in Y 7 for geodesies of bounded length leading to Yp. So we 
have an intersection of set of codimension-one and a set of dimension zero 
dictating the existence of a geometric geodesic through these three points. 

The following result (whose proof is contained in Appendix [A]) justifies 
our use of the term "conjugate" above. We let V&(X) denote the space of 
"b-vector fields," meaning those tangent to dX. 

Proposition 6. The coisotropic manifolds J 7 ^ andJ-t intersect non-transversely 
if and only if there exist cone points Y a ,Yp, a geodesic j(t) with 7(0) e Y a , 
7(t') 6 Yg (with t' < s + t) and a normal Jacobi field W along 7 with 
WeV b (X). 

Remark 4. 

• This result is equivalent to the same statement for the manifolds 
x J 7 a , x J 7 a , as the incoming and outgoing components of J 7 ^ are each 
naturally identified with M.t x ~F^. 

• The transverse intersection stated here is equivalent to the trans- 
verse intersection of the isotropic fibers of the coisotropic manifolds, 
as will be seen in the proof below. In particular, a point of intersec- 
tion automatically entails that the tangent spaces of both manifolds 
contain the tangent to a geodesic and the radial vector field; in both 
manifolds, the space spanned by these two directions is the base of 
the coisotropic foliation. 

• The usual description of conjugate points involves the vanishing of 
the Jacobi field at the endpoints of the geodesic; here, by contrast, 
since W e Vb(X), we find that the metric length of W vanishes at 
the endpoints, hence the variation along W should be regarded as 
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an admissible one for a one-parameter family of geodesies: varying 
the endpoint in Y a should be regarded as varying the direction of 
departure from the cone point, with the location of the end fixed at 
the "point" Y a . 

Finally, we note that it always suffices to prove the results of this paper 
with x replaced by \ which is 1 on suppx; also, by the convexity of the 
Euclidean ball, the hypotheses of the theorems are still satisfied if we replace 
our given K by a larger compact set. Therefore, we will assume without 
loss of generality that our manifold is Euclidean on supp(l — \) an d that 
supp x c K. We also use fi to denote a small neighborhood of the compact 
set K. 

3. Decomposition of the wave propagator 

Let 

L = mm d(Y a ,Yi3) 
denote the minimum distance between cone points. Fix 

Sa « 1, 6* < A 

We let V« be cutoff functions, each equal to 1 in a small neighborhood 
x < 8^/4 of a single cone point Y a and supported in x < 5^; let T e C°°(X) 
equal 1 outside Jl and vanish on suppx; and finally let Aj (j = 1, . . . , N) 
be a pseudodifferential partition of I — ^ip a — Y in which each element 
is microsupported on a set of diameter less than 5a (with respect to our 
arbitrary but fixed metric on S*X) and in particular is microsupported over 
K. We can then further arrange that 

N 

j = l a 

(The error term in fact can be taken to lie in ^ , ~ 00 (X°), which is to say it is 
a smoothing operator with Schwartz kernel compactly supported in X°, and 
thus in particular vanishing in a neighborhood of all boundary components.) 

We will consider sequences of wave propagators sandwiched between var- 
ious of the Aj's: with J = (jo, . . . , jfc+i), set 

Tj = A^U^AjMh) ■ • • A jk U(t k )A jk+1 . 

Associated to each such propagator is a word joji ■ ■ ■ jk+i referring to a 
sequence of elements of the cover, as well as the additional data of a time 
associate to each pair of successive letters jeje+i- We say that such a word 
is geometrically realizable ("GR") if for each pair of successive sets jeje+i, 
in the word, there exist pi e WF' Aj e and e WF' Aj e+1 with 

G,t( 

Pe+i ~ Pi- 
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We call the word diffractively realizable ("DR") if instead there are points 
with 

D.ti 

Pe+i ~ Pi- 

Note that any GR word is DR, and also that the words are read from right to 
left in order to conform with the composition of operators. We also describe 
individual pairs of successive letters in a word as diffractive or geometric 
interactions depending on whether these two-letter subwords are GR or 
DR. We say that such a successive pair of letters interacts with a cone point 
if there exists a diffractive geodesic of length ti between pi e WF' Aj t and 
pi + i e WF' Aj e+1 that passes through some boundary component Y a 

The following result comes direction from ("diffractive") propagation of 
singularities, Theorem 1.2 of [23J; its content is simply that singularities of 
solutions to the wave equation propagate along diffractive geodesies: 

Proposition 7. If the word J is not DR then Tj e TZ. 

Lemma 3. If 5a is chosen small enough and to + • • • + t& > 2To, the word 
J — (jo • • • jk+i) cannot be GR. 

Proof. Shrinking 5a = 2~ m J 0, if there are GR words with fixed time 
intervals, in the limit there must be a geometric geodesic of length 2Tq that 
remains in fi, which is ruled out by assumption. □ 

We also have the following, more granular, reformulation of Assumption [2j 
which we will need in what follows: 

Lemma 4. If 5a is sufficiently small, there do not exists words of the form 
ijk£, where jk is GR and ij, jk, and ki all interact with cone points. 

Proof. If the result fails, then taking 5a = 2~ m gives a family of broken 
geodesies given by the concatenation of 7™, 7™, 7™ where 

• 7o* starts and 7™ ends at a cone point, 

• The end of 7™ and the start of 7™ 1 are within distance 2~ m of each 
other in S*X°, 

• 7™ undergoes geometric interaction with a cone point. 

Then taking m — > 00 would yield (by compactness of dX and Proposition [3]) 
a limiting geodesic passing through three cone points, interacting geometri- 
cally with the one in the middle; this contradicts Assumption [2j □ 

The crucial ingredient in the proof of Theorem [2] is the following lemma 
which shows that the propagator eventually locally smooths data microsup- 
ported in WF'A™ by (n — l)/2 — derivatives. We cannot quite just add 
together all these terms for differing m and conclude the theorem, as that 
still leaves singularities starting in supp^ a to be dealt with; however, we 
will see later on that the singularities near the cone points in these latter 
terms can be moved away from the cone points by applying the propagator 
for short time. 
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In order to make iterative arguments, it is convenient to keep track explic- 
itly of singularities that leave never to return: we let O (for "outgoing" ) 
denote the subset 

= {(t,z,T,()eT* (M\0) : (z, C/r> > 0} . 

(Here we have abused notation by identifying X\Cl with a subset of Eu- 
clidean space to which it is isometric.) The set O is mapped to itself by 
positive time geodesic flow; moreover, any bicharacteristic starting in supp \ 
that escapes Q lies in O over X\Q. We let L 2 H s (0) denote the space of 
distributions that are microsupported in O (hence in Poo on K) and lie in 
I? H s (where as before we use the notation L 2 H s for L 2 ([0,T]; H s )). Let 
L 2 8 S {0) = 1? {H s {0) © H*-\0)). 

Lemma 5. There exist 5a, 5^ sufficiently small that for each m and s, for 
all t > 5Tq, 

U(t)A m : £ s — > L 2 f s+ ( n _ 1 )/ 2 -o + L 2 £ s (0). 

Since distributions in L 2 £ s {0) are smooth on Q this implies in particular 
that xU(t)A m : £ s -> £ s+ ( n -i)/2-o- 

Proof. If all diffractive bicharacteristics starting in A m escape to O in time 
less than 5 To, then the result holds, by propagation of singularities. If not, 
some bicharacteristic must hit a cone point within time To — none can do so 
in longer time, as otherwise it would be a (trivially) geometric bicharacter- 
istic remaining in Q for time To, contradicting Assumption [T] We now let so 
be the time at which the first cone point is reached under the flow. In time 
so + 35i/j, then this particular bicharacteristic is at least distance 25^ from 
the boundary; hence if 6a is small enough then taking to = so + 36^ we find 
that any singularity starting within distance 6a of this one is propagated by 
U(to) to lie at distance greater than 5^ from the boundary, by propagation 
of singularities. (Note that we choose 6a small enough that bicharacteristics 
of length close to to starting from WF' A m interact with only a single cone 
point.) Thus, either U(5To)A m has range in £ s (0) or there exists to < Tq 
such that: 

U{t )A m = Y i A e U{t )A m mod K. 

t 

Now some of the words Im are DR in time to and others (most) are not. 
Those that are not give smoothing terms, so we discard them. For those 
that are, we repeat the construction above, twice, starting in WF ^ and 
WF Ak successively instead of in WF' A m . We may thus write U (5To)^4 m as 
a sum of terms 

U{hTo)A m = ^ U(5T -to-h- t 2 )A J U(t 2 )A k U(t 1 )A l U(t )A m + E + R 
where all words jkim are DR, and where 

E:£ s ^£ s (0), RelZ. 
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The choices of t\ and ty, in the sum depend on k and £ just as our choice of to 
depended on how long it took a bicharacteristic in starting in WF' A m to hit 
a cone point. Since this dependence is not relevant in what follows, however, 
we suppress it in the notation. Note that each tj is less than To, as otherwise 
the bicharacteristics must have escaped to O rather than interacting with 
another cone point. 

Now in each word jkim associated with an element of the sum there are 
three interactions with cone points. Some of these interactions are GR, and 
some are not. We encode this by associating a string of G's and D's to 
a word, so for instance a diffractive interaction followed by two geometric 
gives (reading right to left) the string GGD. We break our sum into pieces 
based on this classification. 

For a word containing two successive -D's, i.e. GDD, DDG, or DDD, 
we claim that the propagator maps £ s — » £ s +(n-i)/2-o- The proof is as fol- 
lows: given initial data in £ s , Proposition [4] tells us that the first diffractive 
interaction results in a solution U{t)A m that is coisotropic with respect to 
J- a (for the relevant cone point a,) relative to H s ~°. This distribution then 
propagates in T*M° so as to preserve this coisotropic regularity (see Propo- 
sition 12.2 of [23] or Lemma 4.7 of [22] from which it also follows). Now 
when the singularities arrive at the next cone point (say, Yg) for the sec- 
ond diffraction, by Assumption [3] (i.e., nonconjugacy) we may apply Propo- 
sition [5] to conclude that the solution is nonfocusing with respect to 
relative to H s+ ^ n ^ 1 ^ 2 ^ . Then the second diffractive interaction puts it in 
£P+(n-i)/2-o^ a secon( j application of Proposition [ZJ 

By contrast a word containing G in the middle, i.e., DGD, GGD, DGG, 
or GGG cannot be realizable by Lemma [4| 

This leaves words of the form GDG as the only remaining summands to 
treat. A geometric geodesic passing starting in WF' A,-, traveling for time 
t2, reaching WF' Ai, and then traveling for time 5To — to — h — t<i must reach 
O since the only other option is for it to reach another cone point, which 
again would contradict Lemma|4j Thus U(5Tq — to — 1\ — t2)AjU(t2)A^ must 
map singularities to O. 

Thus, we have established that every term in the sum representing U (5Tq)A. 
either maps singularities to O or smooths them by (n — l)/2 — deriva- 
tives. □ 



4. Weak non-trapping of singularities 

In this section we prove Theorem [2j which tells us that weak non-trapping 
of singularities holds. 

The proof consists of two steps. To start, we will prove the theorem when 
s = (n — l)/2 — 0. To accomplish this, we apply Lemma [5] as follows. We 
decompose 

(4) X U(t) X = 2 xU(t) X Aj + 2 xU(t) X tpa. 

j a 
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Then for t > 5Tb the first sum has the desired mapping property by Lemma[5] 
since WF' \Aj c WF' Aj. We deal with the second as follows: pick r > 2<L,, 
and smaller than L/50 (recall that L is the minimum distance between cone 
points). By propagation of singularities, if WFu c WF'-^ Q then 

ij)pU{T)v e £qo for all /3, 

hence 

U(t)v -J]^jU{T)v eL 2 £ x + L 2 £ s {0). 

3 

Thus, we may rewrite the second sum in Q as 

2 2 xU(t - t) (AjUWrta) mod K, 

a j 

and again applying Lemma [5j this time with t > r + 5Tb, shows that these 
terms, too, enjoy the desired mapping properties. This concludes the first 
step in the proof. 

To finish the proof, we need to show that further smoothing occurs as 
time evolves. To this end, note that we may iterate the result obtained 
above as follows. Given / = (/o,/i) e £ s we choose xi e equal to 

1 on K and split U(5T ) X f = XiU(5T )xf + (1 - Xi)U(5T Q ) X f so that for 
t > 5T 

U(t)xf = u + v 

where 

u = U(t- 5T )xiU(5T ) X f, v = U(t- 5T )(1 - X i)U(5T ) X f. 

Thus by our previously established results, if xi is chosen with support suffi- 
ciently close to if, u e L 2 ([5T , T]; £ s +{n-i)/2-o) andv e L 2 ([5T ,T];£ s (O)). 
For t > 10Tb we now employ the smoothing result established above with 

u(5T ) = XiU(5T )xf e <? s+ ( n _i)/ 2 - 

now functioning as our initial data (and the previous localizer X replaced by 
Xi) to obtain 

u(t) = U(t - 5T )u(5Tb) 

e L 2 ([10T ,T];f s+2 („_i)/2-o) a+(n-l)/2-0 

Hence overall 

C/(t)x/eT 2 ([10T ,f];£: s+2(n _ 1)/2 _ ) + T 2 ([10To,f];^(0)). 



Further iteration of this argument now yields smoothing by k(n — l)/2 — 
derivatives after time 5/cTq. □ 
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5. Exterior polygonal domains 

In this section we show that the very weak Huygens' principle of Theo- 
rem [2] also holds for the wave equation exterior to a non-trapping polygonal 
obstacle with Dirichlet or Neumann boundary conditions. In particular, we 
suppose that $7 c M 2 is a compact region with piecewise linear boundary. We 
further suppose that the complement M 2 \r2 is connected, that no three ver- 
tices of Q are collinear, and that M 2 \f2 is non-trapping in the sense that the 
doubling described next satisfies Assumption [TJ This assumption is generi- 
cally equivalent to the requirement that all billiard trajectories missing the 
vertices escape to infinity in some uniform time. 

We now form a manifold X by gluing two copies of M 2 \0 along their 
boundaries. This process yields a Euclidean surface (X,g) with conic sin- 



gularities satisfying the assumptions of Section 2.4 (but with two Euclidean 
ends) . As its proof goes through verbatim for a manifold with two Euclidean 
ends rather than one, Theorem [2] then holds for (X, g). 

Suppose now that A is the Dirichlet or Neumann extension of the Lapla- 
cian on M 2 \f2|^] The method of images then shows that Theorem [2] holds for 
M 2 \0. Indeed, by solving the wave equation on the double X and then sum- 
ming (respectively, taking the difference) over the two copies one obtains a 
solution for the wave equation with the Neumann (respectively, Dirichlet) 
extension of the Laplacian on M 2 \f2. 

6. From weak non-trapping to exponential decay 

In this section, we recapitulate the argument of Vainberg [29] as repack- 
aged by Tang-Zworski [28] in the setting of weak non-trapping of singu- 
larities in order to deduce our resolvent estimate Theorem [T] and hence 
exponential energy decay for the wave equation in odd dimensions and the 
resonance wave expansion (Corollary [I]). 

Having established that the weak non-trapping of singularities holds both 
in the settings of manifolds with cone points and of exterior domains to 
polygons, it will behoove us to adopt a formalism for passing from this 
property to resolvent estimates that will simultaneously apply in both cases. 
For this reason we now adopt the "black-box" formalism as used in |28j . 

6.1. Preliminaries. We start by recalling the framework of "black box" 
scattering from Sjostrand-Zworski [26], as used by Tang-Zworski [28] . The 
following presentation follows [28] closely. 

We consider a complex Hilbert space H with an orthogonal decomposition 

n = H Ro @L 2 (R n \B(0,R )), 



2 By an observation of Blair, Ford, Herr, and Marzuola jT], both the Dirichlet and 
Neumann Laplacians are taken to the Friedrichs extension of the Laplacian on the conic 
doubled manifold. 
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where Rq > is fixed. We assume that P is a self-adjoint operator, P : 
H — * H, with domain D c ~H, satisfying the following conditions: 

1r«\Wo) p = H 2 (R n \B(0,R )), 

^-R n \B(0,R )P = A\m."\B(0,Ro) 
(P + i)^ 1 is compact 
PTz-C, C^O. 

Under the above conditions, it is known that the resolvent -R(A) = (P — 
A 2 ) -1 : T~L — > T> meromorphically continues from {A : ImA > 0, A 2 ^ <?(P)} 
to the whole complex plane C when n is odd or to the logarithmic plane A 
when n is even, as an operator from 'Hcomp to T>\ oc with poles of finite rank. 
We denote by T> s the spaces given by (P + i)- s l 2 U. 

The poles (i.e., the resonances of P) wil be denoted Res(P); we count 
them with multiplicity, denoted m(A). We require an additional condition 
on P to guarantee a polynomial bound on the resonance counting function. 

To formulate the additional condition, we use P to construct a self-adjoint 
reference operator P# on 

n* = n Ro @L 2 (m\b(o,r )) 

as in [26j by gluing the "black box" into a large torus instead of Euclidean 
space; here M = (M/iffl)" for some R > Rq. Let N(P, I) denote the number 
of eigenvalues of P& in the interval /. The assumption we need is then 

(6) N(P*, [-C, A]) = 0(A n# / 2 ), A ^ 1 

for some number ^ rajf] 

In addition we have a polynomial bound for the logarithm of the norm of 
the cutoff resolvent in the whole complex plane away from the resonance set 
in odd dimensions and in neighborhoods of the real axis in even dimensions: 
for any X 6 C c °°(M n ), X = 1 near B(0,Rq), 

WxRWxWh^H < C e e c ^\ n * + \ Ae{A:argA<#}\ (J B(Xj,(Xj)~ n *~ e ). 

AjeRes(p) 

Note that the condition Q is satisfied for the polygonal exteriors of Section[5] 
as well as for conic manifolds, with 1~Lr taken to be L 2 (K) in either case, and 
V the domain of the square root of the appropriate Laplace operator (i.e., 
with boundary conditions in the polygonal case and simply the Friedrichs 
extension in the conic case). 



This implies (see |26] and the references of [28]) that the resonance counting function 
N(r) satisfies 

2V(r)= Yi m(\)^Cer n *. 

AeRcs(P) 
|A|<r,arg A<9 
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The wave group of a black box perturbation can be denned abstractly as 
in Proposition 2.1 of Sjostrand-Zworski [27j: 



U(t) = expit Q 



The entries of the matrix representation of U{t) are 

u(t) _ i (DtU(t) U(t) 

where the strongly continuous family of operators U{t) : D s — > T) s+ \ can be 
identified as the solution operator of the following initial value problem: 

{D 2 t - P)U{t)g = for teR 

U{Q)g = 

d t U(0)g = g 

where g e 7i. In other words, 

U(t) = S ^?. 

Since lRi\_B(o,i? Q )^(*)lR n \B(o,Ro) ma P s Schwartz functions to tempered dis- 
tributions we can describe it by its Schwartz kernel U(t,x,y), which is a 
distribution in on R x (R n \B(0,R )) x (R n \B(0, R )). 

6.2. The resolvent estimate. We work in the framework of black-box 
scattering described in Section 6.1 As above, we use the notation U{t) to 
represent the sine wave propagator sin (ts/P^ /yP. In what follows we use 
J 7-1 to denote the inverse Fourier transform: 



^xfW= f e ux f(t)dt 



We use the usual notation in which / denotes J 7-1 / '. 

The main result of this section is an adaptation of an argument of Vain- 
berg [29j as repackaged by Tang-Zworski [28J. It states that the very weak 
Huygens' principle of Theorem [2] implies the resolvent bounds of Theorem [T] 

Proposition 8. Suppose that P is a black-perturbation for which the very 
weak Huygens' principle of Theorem^ holds. Then there exists a 5 > such 
that the cut-off resolvent 

can be analytically continued from Im A > to the region 

ImA > -5 log Re A, ReA>(5 _1 
and for some C, T > enjoys the estimate 



IT] ImA| 



in this region. 
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Lemma 6. Suppose that H\ and H2 are Hilbert spaces, and that N(t) : 
H\ — > H2 is a family of bounded operators that have k continuous derivatives 
in t when (eR, depend analytically on t when Re t > T > 0, and are equal 
to zero when t < 0. Suppose that there are constants jo, k ^ jo + 2, and Cj 
so that for all ^ j ^ k, 



dP y ' 



^Cj\t\ j °- j , for Ret > T. 



Then the operator 

N(X) = J^i x N{t) : Hi — H 2 , for ImA>0 

can be continued analytically to the domain — ^ < arg A < % and when 
|A| ^ 1, it satisfies the estimate 

\\N\\ <C j \\\-ie T ^, forj = 0,...,k. 

Proof of Lemma^ This is a straightforward adaptation of the proof due 
to Vainberg |29l Lemma 4 on page 346] (see also Lemma 3.1 of [28J). We 
include it here for completeness. 

The operator N is defined and depends analytically on A when Im A ^ 
and, for ^ j ^ k, 

f 00 , (9 J 
N = (-i\)- J e lXt —N(t)dt. 
Jo oP 

Let r+ denote the contours in the complex t-plane formed by the interval 
[0, T] and the rays [T, T + ice). When j ^ jo + 2, we have 



N 



HA)-'J r e m ^N(t)dt 0<argA<|, 
{-i*)- j f e™^N(t)dt |<argA<vr. 



The first formula allows us to continue N analytically to the half-plane 
Re A > and the second formula to the half- plane Re A < 0. This also 
yields the estimates for j > jo + 2 and therefore for all < j < k. □ 

Proof of Theorem^ We start by fixing x e C^P(X) and s large (s ^ +2 
should suffice) and then fixing To = T s as in the statement of Theorem [2] so 
that for all t > Tq, 

xU(t)x : £r -> £r+s 

for all r. 

In addition to the cutoff function xi = X> we introduce two other spatial 
cutoff functions, \2 and X3> so that Xi e C^(X), P(l — Xi) = ^o(l ~~ Xi)> 
X1X2 = X2, and X2X3 = X3- 
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We introduce a spacetime cutoff function £ so that ( is independent of 
the spatial variables z on if, < ( < 1, and 

C(M) = 

for some Tq ^ To. Finite speed of propagation and our weak non-trapping 
hypothesis imply that 

(1 - QU(t) X :L 2 ^V S 

for all t. 

For any g e I? = T>q, consider QA(t)xi9i which trivially satisfies 

(D 2 - P) CU{t) X g = - (A 2 - P) (l - CMto, 
^(o)x^ = o, 

D t U(0) X g = X9- 

We now define F(t)g = — {pi — P) (1 — QU(t)xg- Our weak non-trapping 
assumption implies that 

F(t)g e C° (R t ; P s _ 2 ) n C^ 2 (R t ; D ) . 

Note that F(t)g vanishes identically for t < Tq and has compact support in 
t for each fixed z (though the size of the support depends on z) . 
We now define -R(A) by 

R{\) = ((H(t)U(t) X ) , 

where H(t) is the Heaviside function. A simple calculation shows that 

Dt [H{t)U{t) X g] = H{t)D 2 t U{t)X9 ~ iS(t) X g, 
and so R(X) satisfies 

(P - A 2 ) R(X)g = IF" 1 ((D 2 - P) CH{t)U{t) X g) 

= IF" 1 (H(t) (D 2 - P) U{t) X g) + X9 
= iF-\F{t)g) (X)+X9, 

where the last equality holds because the support of F(t) is contained in 
t > 0. 

We now write F(t)g = X2F(t)g + (1 — X2)F{t)g and solve an inhomoge- 
neous wave equation on a flat background. In particular, if Ao denotes the 
(flat) Laplacian on R n , we find V(t)g so that it solves 

{D 2 ~A ) V(t)g = (l- X 2)F(t)g 
V(0)g = D t V(0)g = 0. 
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Using the cutoff function X3 (and the fact that xsO- ~ X2) = 0), we observe 
that 



[l-X2)F{t)g 



A Q ) ((l-xa)V(t)g) 
A ) ((1 - X3)V(t)g) . 



{ D t ~ A ) (X3^(t) 5 ) + 
= -[A ,X3]V(t)g+(D? 

We now define i?#(A) as follows: 

i?#(A) = i?(A) + ij£ x ((1 - X3 )V(i)) 
Since P = Ao on the support of 1 — X3, we observe that 
(P - A 2 ) i? # (A)<? = X g + IT' 1 (F{t)g - (D 2 t - P)(l - X3 )^(0ff) 

= X3 + tf^ 1 (^(i)<? - (1 - Xs)(l - X2)i ? (t)s ~ [Ao, X3] V{t)g) 
= X g + iT- 1 ( X 2F(t)g - [A , Xs] V(t)g) 

= X ( I (X2F(i) - [A , Xs] V(t)) } g. 

In other words, 

(7) R*(X) = J R(A)x(l+^ 1 (X2F(t) - [Aq,X3] 

We claim that the term I+iJ 7-1 (x2-P 1 (i) — [Ao,X3] ^ s invertible in a 

logarithmic region, and so the estimates for %i?(A)x will follow from those 
of X-R*(A)x i n the same region. 

The following four estimates, proved below, will justify the above claim: 



[s-2\ 



,[*-2j 

Given these estimates, the theorem holds with R(X) replaced by i?*(A) 
by pl), (10). One may now take 5 < T _1 , and then there is some constant C 



(8) 




Xfi(A) 


L 2 - 




<: 


Cj\X\ J 


-l e T\ ImA| 


j = o, 1 












(9) 




XF{*) 


L 2 - 


^L 2 


<c 


Cj\\\- 


J e T|ImA| 


j = 0, 1, 


(10) 


x(i-X3M«) 


L 2 - 




<; 


Cj\\\> 


-l e T\ ImA^ 


j = 0, 1 






VDj 








(11) 


[A ,X3]V(.) 


L 2 - 


->L 2 


<; 


Cj\x\- 


j e T\lm\\^ 


j = 0, 1, 



so that for ReA > C and ImA > -JlogReA, one has C|A| _1 e T| Im A| < 1/4, 
showing the invertibility of the claimed term in equation ([7]). Shrinking 5 
then finishes the proof. We must thus only prove the four estimates. 

The first two estimates follow from writing out the Fourier transform and 
noting that for z in the support of x, there is some T so that z) = for 
t > T. Estimate (|8j for j = 1 follows directly from the energy estimate, while 
the estimate for j = uses the energy estimate and integration by parts, as 
XjtX = j) te itx The estimate Q follows by the same sort of integration by 
parts argument and the observation that xF(t) is compactly supported in 
time. The lack of smoothness in the t variable prevents the estimate from 
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holding for all j (and is one of the main differences of the set-up here from 
that using in |29j). 

The other two estimates are somewhat more subtle and rely on properties 
of the free wave group. We start by writing 

(12) V(t)g = (1 - X2)(H(t)U(t) X g - H(t)U (t)(l - X2 )X5 + q(t, z), 

where Uo(t) is the free sine propagator, i.e., Uo(t) = sin t ^/Aq/^/Aq . By 
using the equation for V (and that Ao = P on supp(l — X3) and [X)|, £] is 
order 1), we see that 



(13) 



-A )g = -[A ,X2] CH(t)U{t)xg, 
q(0,z) = D t q(0,z) = 0. 



The inhomogeneous term in (13) has compact support in both space and 



time and vanishes identically for t < 0. 



We now undertake to show that ( 10 ), ( 11 ) hold by verifying them for each 
term on the right-hand-side of (12), starting with q. If E + is the forward 
fundamental solution for the wave equation on M n , then we may write 

X (z)q(t, z) = -\ x(z)E+ (t - t) * [A , X 2] CU{i) X g eft. 
Jo 

If n is odd, then Huygens' principle implies that XQ vanishes identically for 



large t. Thus estimates (10), (11) hold with V replaced by q in this case. 



If n is even, then Huygens' principle no longer applies but if t is large then 
the inhomogeneous term is disjoint from the singular support of E + and 
therefore for large t, 

f 00 f 1=2 
X (z)q(t,z) = c n x(z)\ ((t - t) 2 + \z - z\ 2 ) 2 ([A ,X2] CU(t) X g(S)) dzdt. 

Jo JR n 

In particular, it is analytic for large t and satisfies the derivative estimates 



(10), (11) for large t as well. 



A similar result holds for the second term in (12), since 

^o(t)(i - X2)xg = E+(t) * (1 - X 2)xg- 



Now to see the estimate (10), first observe that the first term in equa- 
tion (12) does not occur because [Ao,X3](l ~ X2) = 0. Combining the 
estimates on XQ an d Wo(t) for large t with the known C s_1 bounds on V{t) 
we apply Lemma [6j 



The final estimate follows in the same manner, but uses the esti- 

□ 



mate ^8^ to bound the first term from equation (12). 



Appendix A. Proof of Proposition [6] 

In this section we prove Proposition |6j By Remark [4] it suffices to prove 
it for x T a and x P/s rather than T a and Pp. 

In what follows, <f>t : TX — > TX denotes geodesic flow for time t and 7( P) «) 
denotes the geodesic in X such that 7(0) = p and 7'(0) = v. We denote by 
7r the projection TX — > X. 
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The kernel of the pushforward 7r* consists of those vectors in TL^ (TX) 
tangent to the fibers of TX and is referred to as the vertical subspace of 
Tf pv \(TX). The Levi-Civita connection on X defines the connection map 
K : Tfa v \(TX) — > T p X , whose kernel defines a horizontal subspace of 
T(p, v ){TX). The connection map K thus provides an identification of the 
vertical subspace of Tu v \[TX) with T p X. Similarly, the pushforward ir* 
provides an identification of the horizontal subspace with T p X. 

We now recall the following characterization of Jacobi fields on a Rie- 
mannian manifold X from Proposition 1.7 of Eberlein |17j . 

Lemma 7. Ifj(p >v ) * s a geodesic in X then there is a one-to-one correspon- 
dence between Jacobi fields along 7 and vectors in T(~,,v) (TX) . In partic- 
ular, for £ e T(p tV \(TX), let Y^(t) be the unique Jacobi field along 7 with 
^c(O) = 7r *C (the horizontal part of C) and 1^(0) = KQ (the vertical part of 
Q). Then Y^(t) is given by the following: 

In other words, the Jacobi field and its derivative are, taken together, 
invariant under the geodesic flow. 

Note further that because the metric is flow-invariant, if 7r*£ and 
are both orthogonal to v e T p X, then the Jacobi field Y^(t) is everywhere 
orthogonal to n /{t). 

We now abuse notation and consider and T\ as subsets of TX 
rather than T*X. We further abuse notation and use £ and rj to denote co- 
ordinates on the tangent bundle (with ^d x + r]-d y denoting the corresponding 
tangent vector) rather than the cotangent bundle. We note further that they 
extend to the boundary and use a bar to denote their completions, i.e., in a 
neighborhood U of Y a , 

= U n {0 < x < s, yeY a , (el, 77 = 0}. 

Observe that if s > and £ > 0, then S (O, y, £, 0) = (s£, y,£,0). A similar 
statement holds for s, £ < 0. 

We also require the following lemma, which describes how the pushforward 
acts on the tangent space to x J 7 a- 

Lemma 8. Suppose (p, v) = (x, y, £, 0) e lies in a small neighborhood of 
Y a and that £ e T/ p>v -\ { J~a) c T/ pv \(TX). IfK( andir^C, are both orthogonal 
to v then KC, = and 7r*£ = ■ d y . 

Moreover, for small t (i.e., t so that x + t£ ^ and 4>t(p,v) still lies in 
this small neighborhood), then 

K{<j>t)*( = 0, 7T*(^t)*C = C' ■ dy 

The proof of Lemma|8]is a simple calculation and omitted here for brevity. 
We also remark that because X J~^ is given as a flow-out, if £ is tangent to 
this flow-out then so is (cj>t)*C- 

We now turn our attention to the proof of Proposition [6j 
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Proof of Proposition Observe that if x J-& and x J rt p intersect then a seg- 
ment of a geodesic connecting Y a and Yr lies in their intersection, i.e., 
x J 7 a n X F t p contains a segment of a geodesic 7 with 7(0) e Y a and 7(t') £ Yp 
with i' < s + t. 

We start by assuming that x J-& and -^J-^ intersect non-transversely at 
the point (p, v) e TX. We let 7 be the geodesic through (p, v) and observe 
that 7 connects 1^ and Yg. Both flow-outs have dimension n + 1 and the 
tangent space to their intersection contains both the direction of the flow 
and the radial vector field. Because the intersection is non-transverse, it also 
contains a vector £ linearly independent of the previous two. We may thus 
assume that both 7r*£ and KC, are orthogonal to v. By pushing £ forward 
by the flow and applying Lemma [7J we obtain a family of vectors tangent to 
both flow-outs that corresponds to a normal Jacobi field along 7. Lemma [8] 
then implies that this family has a limit at the cone points that projects to 
a b-vector field, since the coefficient of d x vanishes at the boundary. This 
completes one direction of the proof. 

Conversely, suppose that W is a Jacobi field along 7 so that W(0) and 
W(t') are both b-vector fields. By standard arguments we may assume that 
W is normal to the flow. By using Lemma[7j we may associate to TV a vector 
field C along 7 in T{TX). Because it is a b-vector field at the endpoints, 
Lemma [8] implies that it is initially tangent to both flow-outs and so stays 
tangent under the pushforward. Because 7r*£ and KQ are both orthogonal 
to the flow, the tangent space of the intersection along this geodesic must 
have dimension at least 3. Each flow-out has dimension n + 1 and therefore a 
transverse intersection must have dimension 2, implying that the intersection 
is non-transverse. □ 
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